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ABSTRACT 


The  evolution  of  long,  two-dimensional,  internal  waves  of  small 
amplitude  is  studied  for  a  density  stratification  that  excludes  KdV  solitary 
waves.  Experimental  data  are  compared  with  theoretical  solutions  of 
inviscid  and  viscous  models  of  evolution.  The  inviscid  model  is  the  modified 
KdV  equation  for  which  dispersion  and  (cubic)  nonlinear  effects  occur  on  a 
time  scale  even  slower  than  that  of  the  KdV  equation.  The  viscous  model 
is  linear  and  only  accounts  for  damping  of  wave  amplitudes.  It  is  demonstrated 
that  viscosity  dominates  early  wave  evolution  in  the  measured  data.  Ever,  in 
these  experiments,  the  major  source  of  viscous  effects  is  the  interfacial 
shear  layer;  hence,  the  early  dominance  of  viscosity  is  probable  even  for 
geophysical  scale  flows  with  one  of  these  special  stratifications.  We 
also  show  that  the  finite  thickness  of  the  pycnocline  in  the  experiments 
causes  a  significantly  smaller  phase  speed  than  predicted  by  the  theoretical 
models  which  utilize  a  two-layer  approximation.  A  simple  calculation  based  on 
a  model  with  two  homogeneous  layers  separated  by  a  layer  with  a  linear  strati¬ 
fication  accurately  predicts  the  observed  phase  speeds. 
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1.0  Introduction  and  Major  Conclusions 

It  is  well  known  that  long  internal  solitary  waves  may  occur  as 
either  waves  of  elevation  or  depression  depending  on  the  background  distri¬ 
bution  of  density,  (For  example,  see  the  discussion  of  the  Korteweg-de  Vries 
(KdV)  equation  in  Part  I  of  this  series  and  references  cited  there.)  For 
a  density  stratification  separating  these  two  wave  regimes,  the  weak 
quadratic  (nonlinear)  interactions  cancel  exactly  and  no  KdV  solitary 
waves  are  possible.  The  purpose  of  this  paper  is  to  study  the  evolution 
of  long,  two-dimensional,  internal  gravity  waves  of  small  amplitude  in  a  fluid 
with  one  of  these  special  density  distributions .  Theoretical  and  experimental 
results  are  presented  for  a  two-layer  stratification  with  a  small  change  in 
density  across  the  layers.  Here  the  quadratic  interactions  arising  in  the 
upper  and  lower  layers  cancel  exactly  when  the  layer  depths  are  equal.  We  also 
note  that  there  is  a  corresponding  stratification  for  two-layers  with  a  large 
change  in  density  (Djordjevic  and  Redekopp,  1978)  and  with  a  continuously 
stratified  fluid  (Long,  1965  and  Benjamin,  1967).  To  our  knowledge,  this 
is  the  first  experimental  study  of  wave  evolution  in  one  of  these  special 
stratifications . 

Whenever  quadratic  interactions  cancel  exactly,  one  expects  weak  non¬ 
linear  effects  to  be  dominated  by  cubic  interactions.  This  is  indeed  the 
case  for  the  two-layer  system  and  we  show  that  the  modified  Rortewec-de  Vries 
(mKdV)  equation  governs  inviscid  evolution  on  a  time  scale  that  is  evei 
slower  than  that  for  the  KdV  equation.  However,  since  the  inviscid  effects 
of  dispersion  and  nonlinearity  are  now  so  weak,  questions  arise  as  to  whether 
viscous  effects  still  can  be  treated  as  a  small  perturbation  on  the 
inviscid  dynamics,  as  they  are  in  Parts  I  and  II  of  this  series.  In  fact, 
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the  experimental  data  to  be  presented  here  indicate  that  the  damping  effect 
of  viscosity  actually  dominates  inviscid  effects. 

Motivated  by  experimental  data,  we  first  develop  in  §2  a  linear  model 
for  the  slow  evolution  of  long  internal  waves  that  accounts  for  viscous  damping 
but  neglects  both  dispersion  and  nonlinearity.  This  viscous  evolution  model 
is  solved  for  arbitrary  initial  data.  Importantly,  results  here  demonstrate 
that  the  unusual  (dominant)  role  of  viscosity  in  the  experimental  data  is 
not  necessarily  a  consequence  of  laboratory  scale  effects.  The  major  contri¬ 
bution  to  viscous  damping  arises  at  the  interfacial  boundary  layer  rather 
than  at  the  side  walls  of  the  wave  tank.  Hence,  this  effect  persists  even 
in  a  geophysical  setting  where  confining  side  walls  are  absent.  In  §3  the 
mKdV  equation  is  shown  to  be  appropriate  inviscid  model  for  nonlinear  evolution 
of  long  internal  waves  when  layer  depths  are  equal .  Long-time  solutions  of 
the  mKdV  equation  and  its  linear  counterpart  are  presented.  Finally,  experimental 
data  are  presented  in  §4  and  compared  to  solutions  of  both  the  viscous  and 
inviscid  models.  We  demonstrate  that  vi scous  ef fects  dominate  the  early 
evolution  of  the  measured  data  and  that  the  viscous  model  of  §2  provides 
a  good  estimate  of  these  effects.  However,  both  the  viscous  and  inviscid 
(two-layer)  models  fail  to  predict  the  observed  phase  speed  of  these  waves 
accurately.  The  leading  order  predictions  of  these  models,  which  are  identi¬ 
cal,  exceed  measured  speeds  by  18%.  It  appears  that  the  primary  reason  for 
the  slower  observed  speeds  is  the  finite  thickness  of  the  pycnocline  region 
in  the  experiments.  A  simple  calculation  of  phase  speed  based  on  a  three-layer 
model  with  a  linear  stratification  in  the  intermediate  layer  agrees  to  within 
5%  of  the  observed  speed. 
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2.0  The  Viscous  Model 

As  already  noted  in  Part  II,  the  effects  of  viscosity  on  the  evolution 
of  long  internal  waves  in  a  two-layer  system  are  similar  to  those  for  long 
surface  waves.  Primary  viscous  effects  still  arise  at  the  fluid  boundaries; 
however,  now  we  must  include  the  shear  layer  at  the  fluid-fluid  interface. 
Fortunately,  the  interfacial  contribution  is  simple  to  include  when  the  upper 
layer  depth  and  the  lower  layer  depth  h2  are  equal,  i.e.  h^  =  h2  Eh. 

(The  notation  used  here  conforms  to  that  in  Parts  I  and  II.)  To  see  this  first 
consider  a  simple  oscillatory  wave  at  the  interface  with  wave  number  k  and 
frequency  w  =  ck  where 

c  =  (gAh/2)15  (1) 

is  the  linear  phase  speed  for  a  two-layer  model,  g  is  the  gravitational  force 

per  unit  mass  and  A=(p2-Pi)/P2  where  _•  is  the  mass  density.  (We  note  that 

the  effect  of  any  finite  pycnocline  thickness  is  to  reduce  the  phase  speed 

below  that  given  by  (1) .  The  two-layer  model  is  adopted  first  because  of 

its  simplicity  and  will  be  modified  later  to  account  in  part  for  the  finite 

pycnocline  thickness  in  the  experiments.)  We  assume  that  the  wave  Reynolds 

2 

number  satisfies  Re  =  w/vk  >>  1  so  that  viscous  effects  are  weak  and  boundary 
layers  are  thin.  Also,  we  take  the  kinematic  viscosities  v  of  both  the 
upper  and  lower  fluid  to  be  equal  which  is  valid  for  the  experiments.  A 
parameter  which  characterizes  the  magnitude  of  viscous  effects  is 

e1=6/h<<l  5  =  (  2v/io)is  (2a,  b) 

where  6  represents  the  boundary  layer  thickness. 

The  leading  order  effect  of  small  viscosity  is  an  exponential  damping  in 
time  (or  distance)  of  the  inviscid  wave  amplitude  with  a  damping  coefficient 
y(k)  and  e-folding  time  (or  distance)  y  ^(k).  For  long  surface  waves  in 
a  channel  of  breadth  b,  the  damping  coefficient,  say  >,  due  to  boundary 
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layer  dissipation  it  the  two  sidewalls  and  bottom  is  (e.g.,  see  Landau  anc 
Lifshitz,  1959,  page  100) 

Long  internal  waves  experience  similar  viscous  dissipation  from  the  channel 
bottom  and  the  portion  of  the  sidewalls  in  the  lower  fluid  layer.  The  added 
contributions  due  to  the  presence  of  the  upper  fluid  layer  are  found  by 
recalling  from  Part  II  that  the  horizontal  (inviscid)  fluid  velocities  will 
be  equal  in  each  layer  (and  oppositely  directed)  when  the  layer  depths  are 
the  same.  Hence,  the  sidewall  contribution  in  the  upper  layer  is  identical 
to  the  sidewall  effect  of  the  lower  layer.  In  addition,  since  the  inviscid 
velocities  are  equal  and  oppositely  directed  above  and  below  the  interface, 
the  actual  velocity  at  the  (quiescent)  interface  position  is  zero.  The  zero 
velocity  there  is  equivalent  to  a  "no  slip"  boundary  condition,  and  we  can 
imagine  a  (thin)  solid  boundary  at  the  interface  with  identical  boundary 
layers  above  and  below.  Each  of  these  interfacial  boundary  layers  is 
identical  to  the  bottom  boundary  layer;  hence,  the  total  viscous  effect 
for  the  internal  wave  is  found  by  summing  these  different  contributions  which 

leads  to  a  damping  coefficient 

...  _  3b+4h  wl*1  =  3b+4h  vck]*3 

YtK'  2bh  2j  2bh  2  J'  (4) 

(We  have  assumed  a  free  surface  that  is  fully  extensible,  i.e,,  uncontami¬ 
nated,  so  that  no  significant  contribution  to  viscous  damping  arises  there. 
This  assumption  appears  reasonable  for  the  experimental  data  presented  here.) 
We  note  here  that  the  interfacial  contribution  to  viscous  damping  will 
generally  dominate  accounting  for  67%  of  the  tocal  when  b  >>  h  and  exceeds 
50%  for  b  >  4h.  In  the  laboratory  experiments  to  follow  the  interfacial 


contribution  is  57%  of  the  total  based  on  the  above  analysis. 


Now  consider  an  arbitrarily-shaped  one-dimensional  initial  wave 
n(x,t=0)  where  x  is  the  direction  of  propagation  and  t  is  time.  Introducing 
the  nondiraensional  variables 

f  =  </J~  n/h,  x  =  (x-ct)/h,  t  =  ct/6h  (5) 

(for  reasons  that  will  become  apparent  in  section  3)  the  initial  wave 


becomes  f(X/0)  -  fQ(x)  and  can  be  represented  by  its  Fourier  transform 

{CD 

fo<X)e"lKXdX  (6) 


F  (K)  = 
o 


where  K  =  kh.  Assuming  linear  propagation  of  the  initial  wave  and 
neglecting  dispersion  so  that  only  viscous  damping  occurs,  each 
Fourier  component  will  decay  with  a  damping  coefficient  (4)  which  can  be 


written  in  nondimensional  form  as 
1*5 


1 


— 
2hc  * 


(7a, b) 


CJ 


y  (K)  =  6  |K  |  ,  8-  3  (3+4h/b) 

If  F (K,t )  represents  the  Fourier  transform  of  the  wave  system  at  any  time  t, 
then  the  evolution  equation  for  F(K,t)  under  the  above  assumptions  is 


irr-YF  (6) 

which  is  easily  solved. 

The  wave  form  f(X,t)  is  recovered  using  the  inversion  integral.  For 


initial  data  with  F  (-K)  =  F  (K) ,  as  in  the  experiments  the  final  result  is 
_  o  o 


f(X,T)  =  J 


F  (K)  exp (-SK  T )  cos  KX  dK  . 
o 


(9) 


o 

We  note  that  even  though  the  above  analysis  is  presented  in  a  heuristic 
manner,  a  confirming  perturbation  expansion  is  possible  following  Chester 
(1968)  and  Miles  (1976)  who  consider  weak  viscosity,  inviscid  dispersion, 


and  nonlinearity  for  long  surface  waves. 


to  obtain  significant  nonlinear  effects  for  the  case  of  equal-layer  depths. 
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The  necessary  assumptions  are: 


i)  hi  *  h2  ;  h 


ii)  c. 


kh  <<1 


iii)  e,  =  n/h  <<1 


iv)  e  B  e. 


v)  <<  c' 


0(e2) 


vi)  the  wave  is  strictly  one-dimensional 

where  n  is  a  characteristic  wave  amplitude.  With  these  assumptions  and 

invoking  the  Boussinesq  approximation  ( A  -*•  0 ,  g  -*■  <*>  so  that  Ag  is  finite) 

the  leading  order  equations  are  the  same  linear  hyperbolic  system  discussed 

in  52.1  of  Part  I.  Secular  terms  do  not  occur  at  the  next  order  0(e)  but 
2 

do  arise  at  0(e  ),  corresponding  to  evolution  on  the  slow  time  scale 
2 

t2  =  e  t.  The  dimensional  equation  for  the  rightward-running  wave  on  this 


time  scale  is 


1_  _3n  In  (nf  — 

c  3 1  +  3x  3  (h  3x 


.  2  .3 

*  V  Mr  -  o 

6  3x3 


(14) 


To  reduce  (14)  to  standard  form  we  introduce  the  nondimensional  variables 
in  (5),  so  that  (14)  becomes 

f  -  6f2  f  +  f  =0 

T  X  XXX  115) 

Djordjevic  and  Redekopp  (1978)  noted  that  there  are  no  soliton  solutions 
of  (15)  that  decay  for  |  X I  since  the  cubic  nonlinear  term  and  the 

dispersion  term  have  opposite  signs.  The  long-time  solution  of  (15)  was 
given  by  Ablowitz,  Kruskal,  and  Segur  (1979).  For  time  t  sufficiently 
large  the  solution  near  the  wave  front  is  approximated  by 


f(X,t)  ~  (3t 


W 


X/(3t>  +  0 1 (3 i 


(16) 
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where  W(z)  is  the  second  Painleve  function  defined  by  ordinary  differential 
equation 


d2W 

dz2 


zW  +  2W  ; 


the  boundary  condition  for  (17)  is 
W(z)  t  rAi(z)  as  2  +  ■ 


(17) 


where 


r  =  tanh  V. 

In  dimensional  terms  (16)  becomes 


n(x,t)  n.  h 

ct 


x-ct 


h(ct/2h) 


Xl 


(18) 


(19) 


If  the  dimensionless  wave  volume,  V,  is  sufficiently  small,  then  r 
in  (18)  may  be  expanded  in  terms  of  its  argument.  In  addition,  W(z)  reduces 
to  Ai(z)  so  that  (19)  is  approximated  by 


n(x,t)  = 


h¥-  f2h)^ 

»T  lct. 


Ai 


[  x-ct 
l'h{ct/2h)^ 


(20) 


We  note  that  (20)  is  also  the  long-time  solution  of  the  Kdv  model,  which 
reduces  to  a  linear  equation  for  this  special  stratification.  Even  if  the 
wave  volume  is  not  small,  both  (19)  and  (20)  predict  that  the  wave  front 
becomes  self-similar  with  the  same  similarity  variable;  i.e., 


W 


x-ct  j 
h(ct/2h)  3_ 


(21) 


as  t  ■*  "  with  only  the  function  W  differing  between  solutions. 

4.0  Comparison  wuth  Experiments 

In  order  to  test  the  theoretical  results  of  sections  2  and  3,  a  series 
of  experiments  was  conducted  following  the  procedures  discussed  Part  I  with 
h  =  5  cm  and  A  =  0.052.  The  measured  density  stratifications  at  the 
beginning  and  end  of  the  experimental  series  are  shown  in  figure  1.  Note 


9 


that  the  pycnocline  thickness  initially  is  about  lcmbut  increases  to 
approximately  4  cm  before  the  end  of  the  experimental  series.  (The  increase 
in  pycnocline  thickness  results  primarily  from  molecular  diffusion  as  shown 
by  Hammack,  1980.)  The  initial  wave  for  this  series  of  experiments  is 
rectangular  in  shape  with  a  length  of  122  cm  and  ar.  amplitude  of  0.25  cm. 
Temporal  records  of  the  interface  elevation  were  measured  at  five  positions 
downstream  of  the  generation  region  at  3  m  intervals. 

We  now  estimate  time  scales,  based  on  a  dominant  length  scale  of 

-in  —  <v 

kQ  —  122  cm,  and  an  amplitude  scale  of  r;o  —  0.25  cm.  A  characteristic 

(fast)  time  scale  is  -'0_1  =  (k^  c  )  ^  11  secs.  The  three  parameters 

characterizing  viscous  damping,  inviscid  dispersion,  and  nonlinearity, 

for  the  initial  wave  become,  respectively, 

c,  =  ‘  /h  '0.094  t,  =  k  h  -  G.041  e,  =  n  /h  =  0.050.  (22 

i  o'  i  o  J  o 

Hence,  representative  measures  of  the  slow  time  scale  for  significant  viscous 
damping  may  be  taken  as  the  e-folding  time  y  ^  (k  )  =  130  secs  (using  (4) 
with  b  =  39.4  cm  and  v=  0.01  cm2/sec)  or  the  unit  slow  time  scale 
(f ^  =  120  secs.  This  slow  time  corresponds  to  a  propagation  distance 
of  about  10  m  which  is  within  the  15  m  range  of  experimental  measurements. 
Alternatively ,  a  unit  of  the  slow  time  scale  for  nonlinearity  (and  dispersion) 
to  become  important  is  *  1  hour!  Clearly  then,  viscous  effects 

should  dominate  wave  evolution  over  the  range  of  experimental  measurements. 

In  figure  2,  the  waves  measured  at  the  last  three  downstream  positions 
(x  =  9m,  12m,  and  15m)  are  compared  with  theoretical  predictions  based  on 
the  viscous  model  (9) ,  which  is  evaluated  numerically.  The  theoretical 
results  agree  fairly  well  with  the  measured  data  in  terms  of  the  amplitude 
and  length  scale  of  the  dorminant  portion  of  the  evolving  wave.  The  major 
discrepancy  between  theory  and  experiment  occurs  in  the  arrival  times 


of  the  wave  at  the  downstream  positions.  Measured  arrival  times  at  each 


location  correspond  to  an  average  phase  speed  that  is  18*  smaller  than 
that  of  the  theoretical  model;  we  will  return  to  this  discrepancy  shortly. 

We  have  also  plotted  the  measured  wave  data  in  figure  3  in  terms  of  the 
similarity  variables  (13)  for  asymptotic  viscous  decay.  Obviously,  the 
data  have  not  achieved  asymptotic  behavior  even  though  they  do  appear 
to  be  merging  into  a  single  functional  form. 

To  reinforce  our  hypothesis  that  viscous  decay  dominates  both  nonlinearity 
and  dispersion  for  this  experiment,  we  show  in  figure  4  the  measured  data 
at  the  last  three  stations  in  terms  of  the  inviscid  similarity  variables  (21) . 
The  data  exhibit  similar  trends  but  clearly  have  not  collapsed  into  a  single 
functional  form.  Figure  4  also  shows  the  theoretical  solutions  (19)  and  (20 
for  the  asymptotic  structure  of  the  wave  front.  The  differences  between 
the  measured  data  and  the  (similar)  theoretical  solutions  are  striking. 
Predicted  amplitudes  exceed  measurements  by  500%  while  the  predicted  length 
scales  are  about  50%  smaller  than  those  measured.  (As  in  the  viscous  model, 
the  predicted  phase  speed  is  18%  greater  than  that  measured.) 

Both  the  viscous  and  inviscid  models  discussed  above  predict  a  leading 
order  phase  speed  given  by  (1)  which  exceeds  the  observed  speed  by  18%.  The 
disagreement  between  data  and  theory  seems  to  be  primarily  a  consequence 
the  two-layer  approximation,  which  ignores  the  finite  thickness  of  the 
actual  pycnocline  region.  In  fact,  the  experimental  data  presented  here 
were  obtained  within  one  and  one-half  hours  of  the  final  density  stratification 
shown  in  figure  1,  where  the  pycnocline  thickness  is  approximately  4  cm  or 
40%  of  the  total  fluid  depth.  As  a  first  approximation,  the  effect  of  the 


higher  order  effects  of  nonlinearity  and  viscosity. 
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